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We perform an accurate computation of the production rate for primordial black holes (PBHs). 
The reason is that the underlying mass variance had been overestimated systematically, as was 
shown recently. For scale-free powerlaw primordial spectra, and for a Universe with critical density, 
the mass variance is less than 34% of its value thought earlier for the spectral index in the range 
1 < n < 1.3. We then extend our study to spectra with a characteristic scale and find the accurate 
shape of the corresponding mass variance. For a pure step in the primordial spectrum, the step in 
the variance is smoothed around the characteristic scale ks. For a spectrum with large oscillations 
near ks, we find a pronounced bump in the variance. This could yield a significant part of the cold 
dark matter in the form of PBHs with mass M in the range 5 x 10^^ g ^ ^ 10^^ g. 

PACS numbers: 04.62.+V, 98.80.Cq 



I. INTRODUCTION 

Primordial black holes (PBHs) are black holes that result from the collapse of density fluctuations in the early 
Universe This is the only way to produce black holes with mass smaller than about three solar masses (except 
possibly by the production of black holes in accelerators, as has recently been discussed in the speculative framework 
of branes). This makes them interesting for several reasons. First, PBHs with initial mass < A/, 5 x 10^'* g have 
been already evaporated by Hawking radiation or evaporate at the present stage of the Universe. They would play a 
key role in studying quantum gravitational effects of black holes (see, e.g., ||^). Second, PBHs with any mass could be 
an excellent indicator of the conditions in the early Universe. And third, they could in principle contribute a major 
part to dark matter. As no PBHs have yet been observed (or identified as such), existing observational limits can 
strongly constrain scenarios of the early Universe. 

One of the most successful scenarios is inflation (see, e.g., |^ for a review). It leads to a primordial spectrum 
of fluctuations that can serve as the seed for the observed structure in the Universe. Recent observations of the 
anisotropies in the cosmic microwave background (CMB) radiation strongly support inflation. Depending on the 
exact form of the spectrum, the same fluctuations can also lead to the formation of PBHs (see, e.g., |5)). Most 
models of inflation predict a scale-free spectrum with a spectral index n very close to the scale- invariant case n = 1. 
A significant number of PBHs can only be produced for n > 1 (a "blue spectrum" ) , since these values lead to more 
power on small scales (this was first addressed in ) . Observational limits (both from Hawking radiation and the fact 
that PBHs must not overdose the Universe) strongly constrain n. This, therefore, yields a constraint on inflationary 
models that is independent of the cosmological constraints from the CMB and large-scale structure. 

It has been found recently that the production rate of PBHs was systematically overestimated Q . The reason is 
that one has to use the transfer function (the tool to describe the evolution of fluctuations within the horizon) and 
evaluate it at the onset of PBH formation, not at the present time. The qualitative explanation goes as follows: at 
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each stage of the cosmological evohition, the correct use of the transfer function takes into account the dynamics of 
the fluctuations after horizon-entry. This leads effectively to very different spectra on small scales at different times. 
A general formula to remedy this earlier deficiency has been presented in |Q and some quantitative estimates have 
been given. The extension from vanishing cosmological constant discussed in Q to the nonvanishing case has been 
performed in However, an accurate calculation has to be done numerically. This will be done here and will 
constitute the first aim of the present paper. 

Although scale-free models have been successful so far, it is not clear that the true spectrum must be exactly 
scale-free. The inflaton potential may exhibit distinguished features at an energy scale relevant to GUT theories. 
This would lead to models with broken scale invariance (BSI). One simple example would be a jump in the derivative 
of the inflaton potential, interrupting the slow-roll behavior and leading to additional particle creation and formation 
of PBHs. The spectrum of fluctuations resulting from such a jump can be presented analytically in closed form 
It can cope surprisingly well with observational constraints like the possible presence of a bump in the matter power 
spectrum pO| ]. The BSI model contains two additional parameters: the strength of the jump, p, and its location, fcs- 

In wc have already discussed the simple case of a pure step in the mass variance. The BSI model of ||^ is, of 
course, of much more physical relevance. The second aim of our paper is thus to perform an accurate calculation of 
PBH production in this model. In fact, the production rate can be so high that the PBHs could be candidates for 
massive compact halo objects (MACHOs) and thus contribute a significant part of the dark matter in the Universe 

BSI models have the advantage that they are the only realistic models of PBH formation. The reason is that 
scale-free models need a certain fine-tuning to lead to a significant number of PBHs: if the spectral index is only 
a little more than some critical value n*, the PBH production rate would be so high that it were in conflict with 
observation. On the other hand, for n being somewhat smaller than n*, the production rate would be exponentially 
small and one would never have the chance to observe a PBH. The critical value is n* « 1.3 which is slightly bigger 
than current observational limits (see e.g. p^). Even more, a blue spectrum with constant spectral index close to 
n* is inevitably ruled out due to the overproduction of PBHs which have evaporated by now and contribute to the 
diffuse 7-ray background. The situation will be drastically different in a BSI model, as will become clear from the 
results presented here. PBH production from models with a scale was considered before |lj, but we shall 
reconsider it here in the light of our improved formulae from Q]. 

Our paper is organized as follows. In Sec. II we shall briefly review the fundamental formalism to calculate PBH 
formation, see |^ and 0| for more details. In particular, we shall present the crucial relation between primordial 
power spectrum and mass variance. In Sec. HI we then give accurate results for the cases of scale-free spectra, spectra 
with a pure step, and the BSI model of Sec. IV is devoted to the discussion of observational constraints for each 
of these situations. For the BSI model, in particular, these constraints allow the possibility to form PBHs in the 
mass range 5 x 10^^ g ^ ^ 10^^ g- Such PBHs could play the role of MACHOs in the halos of galaxies and could 
contribute significantly to the dark matter in the Universe. 



II. PBH FORMATION 



It was realized some time ago that a spectrum of primordial density fluctuations would inevitably lead to the 
production of primordial black holes. Let us assume for simplicity that a PBH is formed when the density contrast 
averaged over a volume of the (linear) size of the Hubble radius satisfies 6min < S < Smax, where the PBH mass 
MpBH is just the "horizon mass" Mh, the mass contained inside the Hubble volume. Usually one takes 5min ~ \ 
and Smax = 1; however, recent numerical simulations suggest a higher value of 5min ~ 0.7 

Each physical scale \{t) is defined by some wavenumber k and evolves with time according to A(t) = 2Tra{t)/k. For 
a given physical scale, the "horizon" crossing time tk is conventionally defined by fc = a{tk)H{tk). It is the time when 
that scale re-enters the Hubble radius (here we do not mean "horizon" crossing during inflation but after inflation), 
which will inevitably happen after inflation for scales that are larger than the Hubble radius at the end of inflation. It 
is at that time tk that a PBH might form with mass Mnitk)- Clearly, there is a one-to-one correspondence between 
X{tk), Mff{tk), and k. We can of course also take this correspondence at any other initial time ti and relate the 
physical quantities at both times ti and tk- 

Generally, if the primordial fluctuations obey a Gaussian statistics, the probability density p{6), where 6 is the 
density contrast averaged over a sphere of radius R, is given by 



1 f 

p(S) = —= e 2<.^(-R) 
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where the dispersion cr^{R) = ^^^^^ ) is computed using a top-hat window function, 



o-^(^) = TT^ / dkk"^ W^HikR) P{k) 
27r^ Jo 



(2) 



Here, P{k) = (|i5fcp) is the primordial power spectrum. The averages (. . .) refer in principle to quantum expectation 
values, but an effective quantum-to-classical transition takes place so that it is sufficient for our purpose to deal 
with effectively classical averages. 

Usually what is meant by the primordial power spectrum is the power spectrum on superhorizon scales. On these 
scales, the scale dependence of the power spectrum is unaffected by cosmic evolution. On subhorizon scales however, 
this is not the case, and its deformation is correctly described by the adequate transfer function (see Eq. (|^) below). 
Hence, the power spectrum P{k) on sub-horizon scales appearing in (|^) must involve convolution with the transfer 
function. 

The function WTnikR) stands for the Fourier transform of the top-hat window function divided by the probed 
volume Vw = ^ttR^, 

3 

WTH{kR)^ jT-^ {sin kR - kRcoskR) . (3) 

The top-hat window function is the most physical choice to study the formation of PBHs. The reason is that with 
this choice the smoothed density contrast 5 describes directly the average density contrast in the region relevant for 
PBH formation, i.e. a sphere of given radius. The use of a Gaussian window function, for example, would erroneously 
yield too small values for (^^(i?), see for example Fig. 4.1 in the textbook by Liddle and Lyth Q. 

It then follows that the probability (3{Mh) that a region of comoving size R — {H~^ / a)\tf, has an averaged density 
contrast at horizon crossing tk in the range 5min < 5 < Smax , is given by 

o-nitk) 



P(^Mh) = -j= / e d<5 « Jl} e , (4) 

min 

where a^itk) := '^^{R)\t^- 

The expression (^) for (3{Mh) is usually interpreted as giving the probability that a PBH will be formed with a 
mass MpBH > Mnitk), i-e., greater or equal to the mass contained inside the Hubble volume when this scale re-enters 
the Hubble radius. Strictly speaking, however, this is not true, since does not take into account those regions 
that are underdense on a scale Mh, but nevertheless overdense on some larger scale (see e.g. the papers by Carr in 
m for a discussion). In the Press-Schechter formalism this seems to be taken care of in some models by multiplying 
(m) with a factor 2. Fortunately, in most cases (3{Mh) is a very rapidly falling function of mass, so this effect can 
be neglected. In this case, (]{Mh) does give the probability for PBH formation and thus also (at time tk) the mass 
fraction of regions that will evolve into PBHs of mass greater or equal to Mh- 

The main problem in calculating the production rate for PBHs is the correct evaluation of (||), i.e. to find the 
direct relation between the primordial power spectrum and the mass variance. For this we have to evaluate the power 
spectrum at the horizon crossing time tk, giving P{k, tk). Instead of P{k) itself, one also uses 

4(fc,ifc) = ^ P{k,tk) = ^ fc' ^^{k,tk) , (5) 

where $ denotes the gauge-invariant gravitational potential. We will sometimes use the short-hand notation 
5jf{k,tk) = Sjfjtk). Inserting this into (||), one finds the following expression, where all quantities are evaluated 
at the time tk ||, 

a%itk) = a^{k) S%{tk) , (6) 
^ a'{k) k'^\k,tk) (7) 
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with 

a'(fc) = ^H^itk) f " 5l{tkx)T^{kx,tk) W^Hix) dx . (8) 
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The upper integration limit is given by the Hubble radius at the end of inflation, corresponding to the smallest scale 
generated by inflation. The transfer function T{k, t) is defined through 

P{k,t) = ^^^P{k,U)T^ik,t), T{k^O,t)^l. (9) 

Here, U is some initial time when all scales are outside the Hubble radius, k < aH , and we can take U — (the time 
of the end of inflation). 



III. PBHS FROM PRIMORDIAL SPECTRA 
A. General results with scale-free primordial spectra 



In this paper, we will derive accurate results for the production of PBHs from primordial fluctuations spectra 
possessing a characteristic scale. We are thus interested in computing the variance crf^itk) at horizon crossing using 
the power spectrum of the primordial fluctuations P{k) of interest to us. In order to review the formalism and the 
results already obtained earlier, it is convenient to consider first the case of a primordial power spectrum P{k) which 
is scale-free and of the form P{k) — A{t)k^^ on scales larger than the Hubble radius ( "superhorizon" scales). Then, 
from (||,0), we can relate anitk) to the present value Snito) 0, 
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(10) 

(11) 



where Mnitk) denotes the horizon mass at horizon entry of the scale k, and Mniteq) the horizon mass at matter- 
radiation equality. The quantity fco (io) corresponds to the present Hubble radius (the present time). We have tacitly 



assumed in (|10|Jll|) that the PBH forms before the time of equality, tk ^ t^q, and that the universe is first radiation- 
dominated and instantaneously becomes matter-dominated (with A = 0). A more complicated evolution of the scale 



factor ait) will result in different expressions for (hO|jllj). Note that the relationship Mnitk) ex: k~ for tk tgq can 
be made slightly more accurate than that used in ( |10| , |11D [|| . 

From (0) we now find that a is given for a scalc-frec spectrum by the expression 



a'(fc) 



T^ikx,tk) W^ni^) 



(12) 



The quantities Snitk) and fc5$(fc,ifc) scale like oc k^~ , up to the coefficient (^ in our case) which is due to the 
transition between radiation and matter domination. However, due to the non-trivial integral (^2|), this is not the 
case for the quantity unitk)- 

We see from (|6[0) that the function Q:(fc) introduces an additional dependence on k. Therefore, it will have 
two important effects: First, it affects the value of cr'j^itk) when extrapolated from the present values Snito), or 
kg $^(fco,io), back to the formation time. But in addition, it will affect the shape of aj^itk) as a function of Mnitk)- 
For example, as we will see in detail, a pure step in the primordial power spectrum does not translate anymore into a 
pure step in the quantity aj^itk)- This second aspect is particularly important when one considers primordial spectra 
with a characteristic scale, as we do in the present paper. 

The observational input needed in (|lO|) or ( |ll| ) is the numerical value of the quantity Snito), or fcg $^(fco, to), on the 
present Hubble radius scale which is found using the CMB anisotropy data for large angular scales. Once /cq ^^(fco, ^o) 
or equivalently S^iko, to) is a known number, the overall normalization of the spectrum is fixed. We will take ]l^ 

k^ $^(fco, io) = 1-67 X 10"^ X exp[-0.959(n - 1) - 0.169(n - 1)^] , (13) 

which corresponds to a flat critical density universe i^m.o — 1) and we have tacitly assumed here that the primordial 
gravitational waves background is negligible. We have in particular 



*'(feo,io) 
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(14) 
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We should mention here the effect of a cosmological constant, as recent observations strongly suggest that we live 
in a nearly flat universe with J^a.o ~ 0-7 S'Hd ^m,o ~ In such a Universe, the effect of including a cosmological 
constant would merely reduce the mass variance anitk) by about 15% [||. In a more general quintessence model, this 
decrease is mo del- dependent. 



B. Primordial spectrum with a pure step 



We consider now more sophisticated spectra with a characteristic scale for which Eqs. ( |10yil| , |12| ) do not apply 
anymore. For each specific spectrum we have to compute crnitk) numerically, but it is possible to give a rough 
estimate of what will happen for a class of primordial spectra using some simple models. In j7| we have discussed 
the simplest case in which the step occurs directly in the mass variance uuitk)- Here we consider a pure step in the 
primordial power spectrum at k — ks, with the corresponding time of re-entrance satisfying < t^q , 



kl $2(fco,io) 



fcr 



1 for fc < kg 
p^^ for k > kg 



(15) 



where the ratio of power on large scales to that on small scales is given by the parameter p^. Note that we assume 
that aside from the step, the spectrum has a constant spectral index n. 
One finds 
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where the quantity a{k), which is slightly different from the quantity a{k) defined in (^, is given by 
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(16) 
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a^(fc) 



T\kx,tk) W^hIx) dx 



(18) 



T^{kx,tk) W^h{x) da; 



It is seen that the location of the step introduces a dependence on fc. It is also seen that, as must be the case, the 
effect of the step disappears both for kg kg and kg — > 0. The difference between both limits corresponds just to a 
different overall normalization. 



C. BSI spectrum with large oscillations 



We now further generalize the spectra under consideration and consider a primordial spectrum produced in an 
infiationary model with a jump in the first derivative of the inflaton potential V{(j)) at some scale fc^ and corresponding 
value (jjg of the inflaton. The resulting spectrum for such a potential is of a universal form, and an exact analytical 
expression has been derived by Starobinsky ||] . At re-entrance inside the Hubble radius during the radiation dominated 
stage, one has the following primordial spectrum. 



e^^{k,tk) = - F{k) 



tkc ^ tk ^ te 



(19) 
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FIG. 1: The mass variance cr|f(tfe) (dotted line) and the primordial power spectrum k^^'^{k,tk) c>c Sjj{tk) (solid line) for 
the Starobinsky-type BSI spectrum ( |20[ ) are displayed for p = 7.58 x 10~^. The normalization of each curve is arbitrary. 
This particular value of the parameter p corresponds to p^, the absolute minimal value allowed by the constraint set by the 
contribution of evaporated PBHs to the diffuse 7-ray background (see below). 
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(21) 



where $ is as above the (peculiar) gravitational potential and the quantities are the inflaton potential 

derivatives on both sides of the jump. This expression depends (besides the overall normalization) on two parameters 
p and ks- The shape of the spectrum depends only on fc/fc^; the scale ks only determines the location of the step. 
Note that we have 



F(0) 



9Hf 
2Al ' 
9H^ 



2A^ 



(22) 
(23) 



For p > 1, the spectrum has a fiat upper plateau on larger scales, even with a small bump, and a sharp decrease on 
smaller scales, with large oscillations. Hence, compared to a pure step, the spectrum ( ^ ) is a step which is dressed 
up with a rich structure, in particular large oscillations confined to the neighborhood of the step. But one can also 
have the case with p < 1 with an increase on small scales. This will be the case of interest to us, and this possibility 
with the corresponding variance is displayed in Fig. 1. 

The power spectrum at re-entrance inside the Hubble radius is then given by 



PP(A:,ife) = - F(k) 



We have in particular 



for big scales. Using ( [25 
introduce the notation 



A:^ $2(fco,io) = ^S'M = ^ Fiko) ^ ^ F{0) 



(24) 
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^), we can write the normalized spectrum explicitly. For this purpose it is convenient to 
m ^ ^ fik) = 4^ f{k) . (26) 
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One is then finally led to the following result, 

100 



<^Hitk) ^ ^p-^ Sj,{to) f{kx) T\kx,tk)W^H{x) dx 



(27) 



Though we have introduced here the function F{k) in relation to the particular BSI spectrum considered here, it 
actually applies to spectra of arbitrary shape, whether scale-free or not, with F{k) defined through (p^). Indeed, we 
then obtain 



x^ F{kx) T'^{kx,tk) W|^(x) dec 



(28) 



cf. our general formula (0). It is easily checked that the general expression (|2|) reduces to the expressions derived 
earlier when F{k) = Ck'^~^ for a scale-free powerlaw spectrum or when it has a superimposed pure step. Finally, in 
order to study PBH formation, it is physically appealing to express the results of this section in terms of the mass 
Mnitk) instead of the wavenumber fc, as is actually done in (0,|ll[). We have in particular (for tk,tk^ < t^q) 



k 



(29) 



For all the primordial spectra considered in this Section, we will derive now accurate expressions for crfjltk) and then 
evaluate them numerically. 



D. Accurate expression for ajj{tk) 

We finally have to deal with the transfer function T{k',tk) in order to have an accurate evaluation of the integral 
expression for the mass variance <J%{tk)- In the radiation era, Sk behaves like 



sin(^a) 

Ok oc cos(ua) , 

IJ,a 



where the constant fi is given by 
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It is easily checked that (|30| ) has the right behavior on scales larger than the Hubble radius, 

6k oc fc ^ aH . 

Therefore, after matching on large scales, the following expression is obtained, 

9 / sin(csa;) 
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W^Hix) dx , 



(30) 



(31) 



(32) 



(33) 



which is valid for any primordial spectrum, with F{k) defined in ( p9| ) independently of the particular expression (j^ 
for F{k) that we consider here. For a scale-free powerlaw spectrum, F(k) = Ck^^~^, expression (ESh reduces to 



<(4-) = Sfjih) / x" 
Jo 



sin(cs2:) 
9 I cos(Csa;) 



CaX 



W|^(x) dx 



(34) 



By inspection of ( |28| , |33| ), and (|6|, |l2| , |34| ), we have the following identification, 
T\kx,tk) 



9 / sin(csa;) 

~ ( ~ COS[CsX) 
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so that (p3f) and (134) can be written as 



fee 

I'iiitk) = ^ ' ^^(fcx) W^2^(c,x) W^Hi^) dx , (36) 
and 

fee 

f^H(^fc) = '^H(ife) / ' W^^(c,x) VK^ff(a;) dx , (37) 

JQ 

respectively. For a scale- free powerlaw primordial spectrum, the following numerical value for a'^(k) is obtained: 

5.37 < a2(A:) < 6.83 , 1 < n < 1.3 . (38) 

So, even for a pure Harrison-Zel'dovich (n—1) spectrum, there is a decrease on small scales compared to the normal- 
ization used earlier in the literature, with now 

cT%itk) ^ em 6%ito) , n = l, (39) 

because the scaling af^(tk) oc fc"^^ would be inaccurate. The normalization used earlier in the literature corresponds 

Let us remark that in we have taken as a rough estimate T^{x < 1) = 1, T^{x > 1) x^^. We see that this is 
excellent for < a; < 1; however, the amplitude of the oscillations inside the Hubble radius is larger than its value at 
Hubble radius crossing and thus requires the accurate analytic calculation done here. 

We see that actually T'^jkx, tk) does not depend on k. This will be true as long as Sk on all scales relevant in the 
integration is given by (BOl). This is the case for scales of interest here corresponding to PBH formation very deep 
inside the radiation era and Mpbh > M,,, and it would remain true for much lower masses. Therefore, for a powerlaw 
scale-free spectrum, the only dependence on k comes from the upper limit of integration. However, as the integrand 
decreases rapidly, for all PBHs with M > M*, this limit can be safely replaced by infinity. Hence, we conclude that 
for a scale-free powerlaw primordial spectrum, the function a is essentially independent of k. Incidentally, whenever 
a does not depend on fc, our improved normalization is formally equivalent to taking the "old" normalization and a 
new significantly larger Smin- In order words, there is a degeneracy between taking a smaller constant a with fixed 
6min and taking a larger Smin with fixed constant a, provided they satisfy 

(40) 



Smin J old V ^nim / new 

Still, a dependence on k can appear due to a feature in the primordial spectrum, as can be inferred from (p^). But the 
interesting point is that due to the presence of a(fc), the scale dependence of the primordial spectrum is translated in 
a non-trivial way into the scale dependence of the mass variance at horizon crossing. Fig. ^ shows the mass variance 
for the BSI primordial spectrum (p7|), and for the primordial spectrum with a pure step (1^,^^, both with n = 1. 
It is seen that the pure step in the primordial spectrum gets smoothed when it comes to the mass variance. Hence, 
even if one assumes that the PBH mass Mpbh is exactly the Hubble mass Mnitk), a pure step at the scale k = kg 
will boost the production of PBH on a certain range of masses AM in the neighborhood of Mg = Mnitks)- In the 
example shown on Fig. |2| the rise in the mass variance starts already around 0.05 Ms and the asymptotic value is 
reached around 5 Ms- 

As for the BSI primordial spectrum (|2^), we see that the rich structure of the spectrum ( pO| ) has nearly completely 
disappeared due to the effects of the filtering. Nevertheless, and most importantly, what remains is a noticeable peak 
in anitk)- Note that the peak is not found at kg, but at a slightly smaller scale with kp^ak ~ 1-27 kg. In the past, 
this BSI spectrum was used in order to explain the presence of a bump in the matter spectrum on a scale around 
125/i-iMpc 0. Here, we have a completely different possibility: such a spectrum could produce a fairly localized 
bump in the mass variance at much lower scales. 

We will now put some constraints on the parameter values of our spectra. 



IV. OBSERVATIONAL CONSTRAINTS 



There are various limits on the initial mass fraction I3{Mh) of PBHs of mass Mpbh > Mh, most of them being 
related to Hawking radiation. Here, we consider first the gravitational constraint |Q 

PiMH) < 1.57 X 10-1^ {t^) ' (41) 
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FIG. 2: This picture shows the mass variance anitk) for two scale-dependent primordial spectra with the same ratio of power 
between small and large scales, i.e. the same p: for the pure step, Eqs. ( |l6|jl7| ), and for the BSI spectrum, Eqs. (|2^), both with 
n — 1 and p — 0.1. 



which comes from the requirement that the PBH mass density today must not exceed the present (critical) density 
of the universe, ^pbh,o < 1- Such a constraint was derived, for example, in ^. Note, however, that the numerical 
factor in (^ij) makes use of the adiabatic expansion of the universe and all the relevant g factors . We emphasize 
also that in the derivation of (|4l|), the quantity P{Mh) is interpreted as being the fraction of the energy density 
that forms PBH of mass Mh- Actually, we will be interested in several possible values of the quantity flpBH.o- The 
generalized version of (fflj) can be written in the following way 



PBH. 



o(M) h-^ = 6.35 X 10^^ (3{M) 



lO^fg 
M 



(42) 



Since PBHs of initial mass Mpbh ^ will have evaporated by the present day, the tightest constraint is obtained 
for the PBH with smallest mass not yet evaporated, namely Mh ~ 10^^ g. 

In this section, we will derive constraints on the parameters of the power spectra introduced above. In order to 
do this the quantity aiiitk) has to be evaluated numerically, then inserted in in order to find (3{Mh), and finally 
constrained using ( pl] , ^ ). 

In the following we shall use 



0.7 



0.7, 



(43) 



and evaluate (M) at Mh = 10^^ g. 



A. Scale-free spectrum 



The case of a scale- free spectrum was already considered in |Q using an estimate for a{k), namely a{k) w 1.5 for the 
relevant masses Mpbh ~ 10^^ g. The following upper bound on the spectral index n is then found for the parameters 
(H) (Fig. 1 in corresponds to h ^ 0.5, S^in = 1/3): 



< 1.36 



(44) 
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FIG. 3: The quantity (3{M) is shown for the spectrum ( |l6|jl7[ ), i.e. a pure step in the primordial power spectrum itself, and for 
the case of a pure step in aH(tk) (as discussed in [^]). Both correspond to (arbitrarily chosen) parameters n = 1.31, p = 0.6 
and Ms = Mnitk^) = 10^® g. Both spectra shown here satisfy the gravitational constraint, Eq. (^), which applies only for 
M > 10^^ g. 



Performing now a more accurate calculation based on (|3^), corresponding to 

fee 

a^ik) = [ " a;"+2 1^1^ (c,a;) W|^(x) dx , (45) 
Jo 

and for the values given in (^), we refine the upper bound to the value 

n< 1.33. (46) 

We have discussed in 0| the dependence of this result on the parameters given above: the dependence on h is very 
weak, while that on 6min is significant. We should stress that the constraint derived here on scale- free spectra is 
rather academic in the sense that the production of PBHs with mass smaller than M*, whose evaporation would 
contribute to the 7-ray background, would be higher by many orders of magnitude than that allowed. Of course, the 
upper bound we have found using the gravitational constraint alone could be relevant if the mass scale of inflation 
corresponds to M*. 



B. Spectrum with a pure step 



In the following, we shall restrict ourselves to the case p < 1, i.e., more power on small scales. This is the physically 
interesting case for us, since it offers the possibility to produce more instead of less PBHs on small scales, which can 
have interesting cosmological constraints. Furthermore, the opposite case leads to interpretational problems for the 
quantity /3(M): For p > 1 it is no longer a monotonically falling function of mass and therefore cannot be interpreted 
anymore as the fractional mass density of PBHs of mass larger than M. For such a case, one must be cautious with 
the interpretation of the Press-Schechter formalism. 

Fig. 1^ shows explicitly the important consequence of the non-trivial relationship between the quantities anitk) 
and Sniik) mentioned above. Even though the primordial spectrum has a pure step at some scale k = kg, the 
resulting quantity anitk) does not. Indeed, the convolution of (5/f(ife) with the scale dependent function a(k) leads to 
a "smoothed-out" step. As a result, this "smoothing" will also apply to /3(M), the variable depicted in Fig. ^. This is 
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FIG. 4: Dependence of Umax on p for a spectrum with a pure step in anitk) or in Snitk) at Ms = Mniiks) ^ 10^^ g- For p = 1 
the result for the scale-free case is recovered, whereas for p < 1 the constraint on n is considerably strengthened. 



important whenever one is willing to overshoot the PBH production at some given scale: It places an avoidable limit 
to the sharpness of the "step" , and clearly of any feature, that will survive in aH{tk)- Still, as can be seen from Fig. ||, 
the "step" remains fairly well localized in the neighborhood of Mg (corresponding to 10^^ g in the case displayed on 
Fig. I). 

A pure step in UH{tk) itself has been considered before and leads to a relationship between p and the maximal 
allowed value for the spectral index rimax as depicted in Fig. ^. Note that the result is independent of the scale (as 
long as Ms > 10^*^ g) and that for p = 1 the result of the scale-free case, Umax ~ 1.33, is recovered. No information 
about the allowed height of the step can be obtained if Mg < 10^^ g from the condition (|4|) for a primordial spectrum 
for which the quantity (Jnitk) has a step structure, since by definition the latter constraint applies only for Mg > 10^^ 
g, which is the condition that non-evaporated PBHs are affected by the step. As becomes clear from Fig. ^ for 
Ms ^ 10^^ g one obtains the same results for both spectra, i. e. the relationship between rimax and p as depicted in 
Fig. is independent of whether one has a pure step in an or in 6h itself. In the latter case, however, one also finds 
a (weak) constraint on the allowed height of the step for 5 x lO^'^ g < Ms < 10^*^ g. This is due to the "smoothing 
out" of the step. 



C. BSI spectrum 

Computing (3(Mh) for the BSI spectrum (^), one discovers that it is no longer a decreasing function of mass, but 
that it has a peak at Mpeafc instead. This is illustrated in Fig. || and constitutes the most interesting result of our 
BSI spectrum in the context of PBH formation. As can be seen from the comparison with a primordial spectrum 
possessing a pure step structure shown in Fig. ||, the presence of a bump in aH{tk) is a consequence of the oscillatory 
behavior of the BSI spectrum in the vicinity of ks ■ Let us note that there is a one-to-one correspondence between 
kpeak and ks for each p, the ratio is approximately given by fcpeafe ~ 1-27 fc^, corresponding to Mpeak ~ 0.62 Ms, for 
the values of p of interest in the following, namely p < 0(10"^). 

A most interesting application of the bump is the possibility to produce a significant fraction of the cold dark 
matter (CDM) in our Universe [0. We shall use the gravitational constraint corresponding to flpBH.o < 0.3. This 
corresponds to a conservative upper bound where in our scenario all of CDM is made up of PBHs. The condition 
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FIG. 5: The quantity I3{Mh) is shown for the BSI-spectrum ( P7|) containing a jump in the infiaton potential derivative for the 
parameters Mnitks) = 10^* g, and from bottom to top p = 10 , 9 x 10~*, 8 x 10~^. As can be seen, P{Mh) acquires a well 
localized bump in the vicinity of Ms = Mnitk^)- Note that we have (3- = /3{Mh <S Mpeak) ^ lO"^'^ x (3{Mpeak), an inequality 
sufficient to avoid the severe constraint from the contribution to the 7-ray background of evaporated PBHs with Adn < M,. 



^PBH,Q < 0.3 then implies for (3{Mh)'- 

PiMH) < 0.48 X 10-1^ {-^j ■ (47) 

The strongest gravitational constraint is obviously obtained at Af = Mp^ak due to the well localized bump. Further- 
more, essentially all of the PBHs are produced with M Ki Mpeak- Hence our gravitational constraint has to be written 
as follows, 

^PBH.O ~ ^PBH,o{Mpeak) < 0.3 . (48) 

As can be seen from Fig. |[ for a given Ms, taking smaller values of p increases P{Mh) and therefore also flpBH,o- Thus 
for each Ms there exists a minimum value for p, Pmin, in order to achieve the gravitational constraint ^pbha ~ 0.3. 
In this way we obtain a relation between Pmin and Ms as displayed in Fig. ^ by the solid line. A shift in the location 
of the characteristic mass AIs (and in the location of the corresponding scale ks) to higher masses Ms (equivalently 
lower ks) still much smaller than IO^^Mq is allowed, but lower values Pmin are then needed in order to achieve 
^PBHfl — 0.3. 

It is possible to use the observed extragalactic diffuse 7-ray background to achieve an interesting window for the 
characteristic mass scale Ms of our model. For this we are considering the PBHs which have evaporated by now and 
may contribute to the observed 7-ray background. Their mass lies in the range AM, with 2 x 10"'^^ g < Mh < 5 x 10^'' g 
p9| . The 7-ray background constraint comes from the requirement that the 7-ray flux from PBHs whose mass is in 
the range AM cannot be larger than the observed diffuse 7-ray background, leading to Qpbh,o{^M) < 10~*. It is 
for this reason that the sudden and substantial drop of P{Mh) on smaller scales Mh <C Mpeak is crucial. From (|4^), 
it is easy to obtain the following equality, 

^PBH,o{Mh < Mpeak) P- ( Mpeak 



(49) 

^PBHfiiMpeak) Ppeak \ Mh J 

with /3_ = I3{Mh < Mpeak) and [Speak = P{Mpeak). 

Taking ^pBH,o{Mpeak) — 0.3 and using the 7-ray constraint ^pbh,o{Mh ^ Mpeak) < 10^^, ( ^ ) leads to 

3^3X10- (#^1' (50, 



(^peak \ -^3 



peak 
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FIG. 6: The allowed region in parameter space {p,Ms = Mnitk^)) is shown. The solid line represents those points for which 
^PBH,o{Mpeak) = 0.3. Bclow the solid line we have O.PBH,o{Mp^ak) > 0.3. In other words, the gravitational constraint at 
Mpeak is violated and this region is therefore excluded. The two lines parallel to the solid line represent, from bottom to top, 
those points for which QpBH,o{Mpeak) ~ 0.1 and 0.01, respectively. Below the dotted line, the 7-ray background constraint is 
violated. It is seen that for Ms > lO'^^g, the allowed parameter values p yield a density ilpBH,o{Mpeak) < 0.3 which becomes 
rapidly negligible with growing Ms- 



where Mh is any mass lying in the range AM 

Taking smaller values of p increases the ratio „ 

Pi 



f3- 



as can be seen on Fig. 0. Indeed, using (||) we obtain 



Ppeak 



■ exp 



52 . 

mm 



1 - 



a 



(51) 



peak 



where cr_ = crnitk <^ ^k^^ak) ^peak = <^H(tk^^ak)- Though — — is an increasing function of - we only consider 
here p < 1 and larger p corresponds to a smaller step - the decrease of tT_ itself with increasing p produces a faster 
decrease of the exponent, and hence a decrease of the right-hand side of (|l|). In other words, a smaller p, which 



corresponds to a larger step, increases the quantity 



On the other hand, as can be seen from (pO[) 



"7^ . v^A-L ±±c*j±Avi, cjjo iiL'wi/, "7^ cannot 

Ppcak ^1 1 Pptak 

grow arbitrarily. Hence, this results in an absolute lower bound p-y for the values of the parameter pmin for which 
^PBH,(){Mpeak) — 0.3 is Satisfied, namely p^ ~ 7.58 x 10^'', and a corresponding upper bound on Alg at 5 x lO^^g. 
For higher Ms, Pmm < P-y is needed in order to achieve Qpbh.o = 0-3i E^nd therefore the 7-ray background constraint 
is no longer satisfied. All these constraints are displayed in Fig. |[ Of course we obviously obtain a lower bound on 
Ms around 4 x lO^^g (for the same gravitational constraint) because we would violate again the 7-ray constraint due 
to the presence of the bump close to the mass range AM. To summarize, requiring that ^lpBH,o = 0.3, the 7-ray 
background constraint can be satisfied for a range of characteristic masses M^ with 



lOM, < Ms < 10^' g 



(52) 



An observational constraint on the allowed range of MACHO masses could therefore give a severe constraint on 
our model as it must correspond to the window (p2|). Though such observations do not exist for the moment, it is 
interesting that our model has this built-in constraint. In particular, it is interesting that our model can be easily 
discriminated from models where PBHs with M ^ Mq are produced. It is instructive to study the dependence of the 



window ( |52D on the present density of PBHs. We can relax the gravitational constraint yielding flpBH.o = 0.3 in order 
to allow a bigger upper bound on Ms and to enlarge the window ( |5^ ) . Actually, we can consider the 7-ray constraint 
independently of the present density of PBHs, i.e., independently of the value (3peak- Obviously this constraint is 
given by 



14 



I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 



I I I I I I I I I I I I I I I I I I I 




p X 10" 



FIG. 7: The ratio 



P- 



with Ppeak = l3{Mpeak), is shown as a function of p. For p > 7.5 x 10 " the ratio is sufficiently 
small in order to satisfy the severe upper bound set on the contribution of evaporated PBHs to the extragalactic diffuse 7-ray 
background. 



with Mh G am. The constraint (|5^) results in a value for which is independent of Alg as long as Mg > 2.5 x 10^^ g. 
Hence for all those Ms values, the minimum value in order to satisfy the 7-ray constraint will be the same, namely 
the value mentioned earlier, p.y = 7.58 x 10^*. Thus, in the way shown in Fig. ^, relaxing the requirement on the 
present density of PBHs, for example setting flpBH.o = 0.1, allows a slightly weaker constraint for the upper bound 
on Ms ~ 4 X 10^"'^ g, and one gets a maximum Mg ~ 4 x 10^'^ g by requiring only ilpBH,Q = 0.01. The lower bound 
on Ms remains essentially the same around 3 x 10^^ g. 

We would like to make here a final comment. Recent numerical simulations seem to imply that PBH formation 
proceeds via near-critical collapse [|l^ whereby PBH of different masses around the horizon mass Mnitk) = Mp^ak 
could be formed at the same time tk, with 



M = K[5 - 5cV 



(54) 



For 7 = 0.35, following l20|, the following result is obtained. 



PBH 



{M,tk,. 



dM 



This is turn yields 



3.86 



peak 



M 



peak 



2.86 



exp[— 1.35 



M 



peak 



2.86 



(55) 



PBH,tot 



(tfc^^^J = 0.8 /3(Mpeafc) 



(56) 



a result similar to that obtained in the more traditional view, where PBHs with different masses M are produced at 
different times tk with M — Mnitk)', in the latter case one obtains ^pBH.tot{tkpcak) — Ppeak- Hence, our constraint on 
^PBH,o{Mpeak) using (^ is stronger than the requirement ^lpBH,totito) < 0.3 based on (p6|). Note that Mp^ak {Mh 
in the notation of ||20| ) is the Hubble mass corresponding to the peak in the quantity /3{M), not in the primordial 
spectrum itself, because one has to distinguish between the quantities (JH{ik) and Snitk)- The above mentioned 
approach does not take into account the properties of the feature in the primordial spectrum, in particular the width 
of the corresponding bump in (3, besides the value (3{Mpeak) which serves as a kind of overall normalization of the 
PBH abundance. On the other hand, it does account in a consistent way for the production of PBHs also in the 
monotonically increasing part of P{M) around Mp^ak, M < Mp^ak- Clearly, the more spiky the primordial spectrum, 
the better this approach is expected to be, so our model provides such a concrete spectrum based on an underlying 
inflationary dynamics. An accurate determination of the mass spectrum would therefore yield valuable information, 
not only on the underlying primordial perturbations spectrum but also on PBH formation itself. 
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V. CONCLUSIONS 

New expressions were derived recently for the improved computation of the mass variance at horizon crossing 
. We had found that for a scale-free powerlaw primordial spectrum, the mass variance for PBH masses around 
Mh, « 5 X 10^^ g is significantly lower than it had been assumed in the literature. In this paper we have refined 
our estimate and confirmed this discrepancy: We find that 5.37 < a^(M*) < 6.83, for 1 < n < 1.3, see Eq. (|3^), 
in particular <J^{ti.) — 6.63 (5|j(to) for & scale-invariant (Harrison-Zel'dovich) primordial spectrum. For scale-free 
primordial spectra, the quantity a is further independent of M for a wide range of scales. This is to be compared 
with the value cr|^(ife) « 25 5jj{to) which corresponds to the normalization adopted previously in the literature. 

We have further applied our formalism to primordial spectra which are not scale-free. Two specific spectra possessing 
a characteristic scale were considered. The first one is purely phenomenological and corresponds to a pure step in 
the primordial spectrum. We find that the step is translated into the probability P{M) in a non-trivial way: It 
acquires a smoothed step centered around the characteristic mass Alg. The second primordial spectrum corresponds 
to a BSI spectrum which is obtained from inflationary perturbations when the inflaton potential has a jump in its 
first derivative at some scale ks ■ It would of course be of interest to relate this scale to a scale of fundamental physics 
at the GUT scale. In this case, we find that the mass variance cr'jj{tk) displays a sharp, pronounced peak centered 
at the scale kpeak ~ 1-27 kg which is not the exact location of the maximal peak in the primordial spectrum itself. 
Again this is due to the non-trivial relationship between the primordial perturbations spectrum itself and the mass 
variance at Hubble radius crossing a^{tk). This in turn yields a pronounced bump in the quantity (3{M) around 
Mpeak which has interesting applications. We have considered here in detail one such application ||ll[, the possibility 
to produce in this model a significant part of the cold dark matter in our Universe. Interestingly, we find that the 
mass of the produced PBHs is constrained to be much smaller than one solar mass. Indeed, the characteristic mass 
of the spectrum, essentially the same as that of the produced PBHs, cannot be fixed arbitrarily because of the severe 
constraint coming from the 7-ray background. The inclusion of a non-negligible spectrum of primordial gravitational 
waves and running scalar and tensorial spectral indices and ut (see e.g.[^l|) could enlarge the class of models 
that can realize our scenario. A thorough study should probably await a more careful determination of Smm and, 
to less extent, of the cosmological parameter J^a. At the hand of the example studied here, it is clear that spectra 
with characteristic scales on very small scales, much smaller than those probed by the CMB anisotropy or large-scale 
structures, offer interesting new avenues. 

We have studied here concrete models with concrete assumptions, for which we have obtained precise quantitative 
results. It is clear that these results may be different for other models. They may also depend on the assumptions 
made. In particular, the infiuence of non-Gaussian fluctuations could strongly change the PBH formation rate p5| . 
However, all situations can in principle be dealt with using the formalism presented here. 
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